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ABSTRACT 



The general solution of the behavior of a viscously damped vibrating 
system having several degrees of freedom, as developed In the literature 
via matrix methods. Is treated In detail here and made the basis of a 
digital computer program which Is capable of determining the natural 
frequencies, the mode shapes, and the displacements of each mass as 
functions of time. This program, which Is written In FORTRAN 60 language 
for the Control Data Corporation 1604 computer, affords several output 
options. It does not treat cases of supercritical damping or cases In 
which two or more natural frequencies are the same. 
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Nomenclature 



A independent variable in reduced system 

A* real part of a complex matrix 

B independent variable in reduced system 

B* imaginary part of a complex matrix 

C damping matrix 

D driving force amplitude 

DI imaginary component of driving force amplitude 

DR real component of driving force amplitude 

F(t) forcing function 

G column of constant coefficients 

I identity matrix 

K stiffness matrix 

M mass or inertia matrix 

N degrees of freedom 

0 null matrix 

p eigenvalue 

R independent variable in reduced system 

r,s subscripts designating different modes 

SS steady state solution 

TS transient solution 

U eigenvector 

V real component of eigenvector 

W Imaginary component of eigenvector 

X generalized displacement, a column vector 

• 

X generalized velocity 
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generalized acceleration 
dependent variable in reduced system 
dependent variable in reduced system 
decrement factor ( So>n) 
damped natural frequency 
dependent variable in reduced system 
excitation frequency 
natural frequency of rth mode 



CHAPTER I 



INTRODUCTION 

1.1 General Remarks . The analysis of a subcritically damped multiple 
degree of freedom vibrating system necessitates obtaining the solution of 
a complex eigenvalue problem to determine the natural frequencies and mode 
shapes of the system. Although the analysis presented in the literature 
for systems with two degrees of freedom may be extended to systems with 
more than two degrees of freedom, manual calculations are too laborious 
to be practical. Therefore the natural frequencies are usually found by 
ignoring the presence of the damping and solving the resulting real eigen- 
value problem. This simplifies the problem considerably and provides a 
good approximation provided the damping is light. Another technique 
employed is to solve the real eigenvalue problem and then account for 
the damping using perturbation techniques. However, even in the absence 
of damping, systems involving more than two degrees of freedom usually 
require iteration or trial and error techniques (such as the Stodola or 
Holzer methods) to obtain the mode shapes. 

The advent of the electronic digital computer has eliminated the 
necessity of ignoring the damping component and increased the size of the 
system for which solutions can be obtained. Although not entirely devoid 
of error the digital computer is highly reliable and its speed of opera- 
tion has made it an invaluable tool in engineering analysis and design. 
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1*2 Method of Analysis * The multiple degree of freedom damped vibrating 
system is described by a set of N linear second order differential equa- 
tions, where N denotes the degrees of freedom involved* Utilizing matrix 
analysis and generalized coordinates the system is described by 

MX + C)< + KX = FOt) 

which is the same form as the single degree of freedom system* However, 

• • • 

M, C, and K now represent square matrices and X, X, X, and F(t) are 
column vectors* The analysis presented in what follows is performed using 
matrix technique because of the compactness of the notation and the order- 
ed computational procedure which is ideally suited for digital computer 
programming* 

The mathematical analysis is demonstrated and equations are derived 
which describe the time behavior of free and forced vibrating systems. 
Finally a digital computer program is presented which performs the opera- 
tions indicated in the mathematical analysis to determine the natural 
frequencies, mode shapes, and time behavior of the damped vibrating system* 
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CHAPTER II 



MATHEMATICAL ANALYSIS 

2.1 Eigenvalues and Eigenvectors . The viscously damped vibrating system 
with one degree of freedom is described by a linear, second order differ- 
ential equation. 

MX +CX t KX^ P(0 

M, C, and K represent the mass, damping and stiffness of the system, and 
F(t) the forcing function. A system with many degrees of freedom is de- 
scribed by a set of second order differential equations similar to the 
case of the single degree of freedom system. Using matrix notation the 
system is described for the case of free vibration by 

(1) MX + C>C t KX = 0 

where M, C, K, are now square matrices of order N, representing, as in 
the single degree case, the mass, damping and stiffness matrices of the 

system. It is always possible to write the equations so that these mat- 

• «• 

rices are symmetric. N is the number of degrees of freedom. X, X, and X 
are column vectors of order N representing displacement, velocity, and 
acceleration in generalized coordinates. 

Premultiply equation (1) by M (Capital letters will hereafter re- 
present matrices and lower case letters scalar constants) to obtain; 

(la) X + H CX + KT'kX = O 

« 

W. J. Duncan and A. R. Collar w have shown a simplified method of find- 
ing the eigenvalues and eigenvectors by writing the second order differen- 
tial equation in reduced form as a first order differential equation. The 

i 

Numbers in brackets refer to references listed in the bibliography 
page 23 
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reduced form is obtained by introducing velocity as a dependent variable. 
Let r -I 



X 



(2a) 




Here Y is a partioned column matrix of order 2N with the N displacement 
components in the upper half and N velocity components in the lower half. 
Let 



(2b) 



A - 



O 1 

-m’'k -m'c 



Here 0 is an nth-order null matrix and I is an nth-order identity matrix. 
The reduced equation then becomes 



(3) Y = AY 



A solution of a first order differential equation such as equation 

(3) may be taken as 

(4) Y* u 

Substitution of the assumed solution into equation (3) gives 

(5) pe^^u = 

Dividing out the exponential term and rearranging, we have 
(5a) 1 a - pi] U = O 

To have a nontrivial solution, it then follows that 
det \k - 

This is the characteristic matrix of the reduced system. The eigen- 
values or roots of the characteristic equation are p^; r = 1, 2, ....2N 
and the eigenvectors are U^; r = 1, 2, ....2N 
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»»here 







The first subscript indicates the position in the column and the second 
subscript corresponds to the r*^^ eigenvalue. 

As long as the damping is less than critical in each mode, the re- 
duced system will yield N complex conjugate pairs of eigenvalues and 
eigenvectors. The eigenvalues will be of the same form as the sub- 
critically damped single degree of freedom system. 

^ + J to,. Vl-Sf 

'pt, - ~Sr<^r - jw,./ 

The bar indicates the complex conjugate and the subscript indicates the 
mode, and where 

- damping ratio 

- natural frequency 

- damped natural frequency 

>^<0,, - decrement factor 

The corresponding eigenvectors are 

Ur ^ Vr-jWr 
2. 2 Homogeneous Equation . 

2.2.1 Orthogonality Relations . In the case of the undamped vibrating 
system it is always possible to choose a set of coordinates in which the 
mass and stiffness matrices are uncoupled. However in the damped vibrat- 
ing system, unless the damping matrix is proportional to either the mass 
or stiffness matrix, a set of coordinates which will uncouple the equations 
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of motion cannot be found without a knowledge of the eigenvalues and eigen- 
vectors of the system. 

K. A. Foss has developed a set of orthogonality relations for the 
eigenvectors of a damped linear dynamic system from which coordinates 
which uncouple the equations of motion may be found. [ 2 ] The technique of 
reducing the second order differential equation to one of the first order 
is again employed, with slightly different notation. 

Let 





‘0 M 


, pM o' 


1 


X 


R = 


_M C 


(6b) B=[^q 


j (6c) 2 = 


X 



Using the above notation, equation (1) becomes 

(7) TR2 BH = O 
Furthermore let 

Prt 

X = 

• ,1 
X = Pr^r'5 

Then 

P ^ 

(8) H = 

PrU; 

Substituting equation (8) into equation (7) and dividing out the exponen- 
tial factor we get, 

(9) Pp R = o 

Since the eigenvalues and eigenvectors occur in complex conjugate pairs 
the complex conjugates of and p^ also satisfy the homogeneous equation. 



where 
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Denoting the complex conjugate by the subscript S, equation (9) may be 
written 



PsR^>s B$s 

Or since M, C, and K are symmetric, R and B must be symmetric, and equa- 
tion (10) may be written 

(10a) P5 §s^R = o 

Premultiply equation (9) by , postmultiply equation (10a) by 

and subtract (10a) from (9). 

- 

Therefore unless r = s the orthogonality relations are 

(11) =o 

(12) ' O 

2.2.2. General solution of the homogeneous equation . In the general 
solution of equation (7), 2N arbitrary constants of integration must be 
evaluated. The 2N initial conditions of displacement and velocity are 
used to evaluate these constants. Assume a solution of the form 

( 13 ) 

where represents the 2N constants to be determined. 
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The initial conditions may be expressed as a vector expansion of the 



natural modes 



(14) Z 



2N 






Premultiply both sides of equation (14) by R 

(15) $5RHW= y 



rs I 



Using the orthogonality relation, equation (11), the summation simplifies 
and equation (15) becomes 



(16) $Irhco) = 



Or 



(l„ 



The general solution is therefore found by substituting the value of the 
constants of integration into equation (13). 



(18) gc°)e' 



I 

It follows from an evaluation of equation (18) at time zero, that 

- 4 Cr_. j 

T = f 

Equation (18) may be written in the following partitioned form using the 
nth-order system notation. 



(19) 



xCt) 







Mr 

_Uv- Ur 



t U^uTnxCo) + PrU^UJCXCol 
r HKCoV -t Pr-UrUZMX6^ U^uJCxfe) J 



*-X \ 



2p,.ujHUi- + C Ur 



from which 



ZN 



(20) XC^V,, 



^ [UrUjtUXCo') 1- pr-UrUjKAXCo) f UrUjCXfijJ ?r^ 



<-st 



1- U7C Ur 
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It should be noted again at this point that for the sub-critically damped 

system being considered, p and U are complex quantities and occur in 

r r 

complex conjugate pairs. The matrix operations indicated in equation (20) 
have been carried out by breaking up the complex quantities into real and 
imaginary parts and then recombining. The operations are relatively 
straightforward but lengthy, therefore only the final form will be pre- 
sented. 

N 

(21) 

rsi 







f fir HrM + S;rC)><^0^ Cos /9, 









f HrC) X 



Me 






Irsl 



Where; 

+• j(Sr Ur«'^r-‘-jWr 

= - Z<Xr C '/r ^ ■ Wr H Wr) - 4/3r V7 M + v7 C V». - wj C Wi^ 

b^» /SrCVr^^Vr -W7MWp)-40^rVr MU4- t ^ VjC Wr 
Ar= VrVp^ -WrW7 
3r = ^rWr 4- WrYr 
Qr- = 0.rAr + bpBr 

Hr ” 



The factor of two in Che above equation results from Che reduction of the 
summation from 2N Co N. When the indicated operations of equation (20) 
are carried out from r = 1 to r = 2N the imaginary components sum Co zero; 
therefore all quantities in equation (21) are real. Equation (21) has 
been given by S. H. Crandall and R. B. McCalley in reference 3. 
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2.3 Non -Homogeneous Equation 

2.3.1 Steady State Solution of the Non-Homogeneous Equation . Once 
again using the reduced form, the non-homogeneous equation becomes 

(22) RH +-BZ. = FCV) 

Where 

- 

In a linear nth-order differential equation the solution of the non- 
homogeneous equation may be found, once the solution of the homogeneous 
equation is known, by expanding Z in a modal series. [ 2 , 4 ] 

Thus , 

(23) = ^ 

I 

here is a variable parameter. 

Upon substitution of equation (23) in equation (22), we obtain 

(24) = FCV) 

t-sl 

-T 

Premultiplication of equation (24) by and making use of the 

orthogonality relations, yields 

(25) 

Or 

(25a) 

Where 

-prRn = 

Fo’ §r^Ck) 



O 

m 
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Using Che Laplace Transform technique^ Che soluClon of equacion ( 25 a) 
is easily deCermined. 

Rearranging 




The Laplace Transform is 

(b) ZVC%) - p^V(s^ ^ ^ 

Therefore 



(c) 



yc%) X 






(5-pr) 

The inverse Cransform of (c) is 



( 26 ) F^(r)d't 

o 

Since = - c<^ + J 

( 26 a) iC^) = VnC't) (coS/Sr (v-Tr) SiM/ 3 r(>-'t))< 5 l't 



SubsCiCuCe ( 26 a) inCo ( 23 ) 

t 

(27) jf (^) (co%f3^ 0 -'t) f J 5 1 M |3 r (V -'t)) d r 

>-:l ^ 

As in Che case of Che homogeneous soluCion Che reduced sysCem quancicies 
are replaced by Cheir nCh-order equivalence. The final form of Che 
sCeady sCaCe soluCion Chen becomes 

( 28 ) cos^Srtt-'e'ia't 

N ^ > 

+ N I fC't) e siM/Sr 

/ Or j 

r=i o 



The notation in equation ( 28 ) Is the same as that used In the homogeneous 
solution. 
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2.3.2 Steady State Solution with Sinusoidal .Excitation . A special, 
but very common problem is the case in which the forcing function, -PC^) » 

is sinusoidal. Expressing the forcing function as 

(29) 

where - indicates **the real part of** 

D - column of driving force amplitudes, possibly complex to admit 
phasing 

CO - excitation frequency 
The non-homogeneous equation therefore becomes 

(30) MX ^ Ck fKX = 

Assume a solution of the form 

(31) x(^') - CR. 

where Q is a coltxmn of undetermined constant coefficients. Substitute 
the assumed solution in equation (30) and divide through by the exponen- 
tial factor to obtain 

(32) -o/MG, +ja.>CG, t KG = D 
Or 

-J 

(33) (S = wC]d 

Therefore the steady state solution for the case of sinusoidal excitation 
is 

(34) X 



-I 






CO I 
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2,4 Transient Plus Steady State Solution . 

2.4.1. Initial Conditions . The total solution is given by 



(35) = T5 t SS 



where TS is the so called transient solution and SS the so called steady 
state solution. However in the total solution the Initial conditions 
of displacement and velocity in the transient solution must account for 
the initial values in the steady state solution. This is accomplished 
as follows. In the case of sinusoidal excitation the steady state 
solution was found to be 



To allow for phase differences in the exciting forces on the masses, let 
D be given by 



(34) 




(36) B = DR + jDI 



and let the inverse of 




be denoted by 



(37) A +j'B 




equation (34) 



becomes 




Carrying out the indicated multiplications, we obtain 
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and the velocity is given by 



(40) X (^) - io (a DR - b'di) s I m to ^ - 60 (b'DR ^ a'di) CoS to ^ 

Therefore the initial values of displacement and velocity of the steady 
state solution are obtained by evaluating equations (38) and (39) at O 

(41) = a'dR-'S^DI 

(42) = - CO (b'dR + a'di) 

The initial conditions used in determining the time behavior of the 
transient portion for the total solution therefore become 

(43) Xe Co) ^ X(g) - (a'-DK - e'Dl) 

(44) V^Co) ~ X(o) 60 (-BDR t A'di) 

The complete time-solution is thus obtained as indicated in equation (35) 
where the transient solution is found by using the modified initial condi- 
tions of equations (43) and (44). 
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CHAPTER III 



PROGRAM DESCRIPTION 

3.1 General Remarks . A digital computer program ''PROGRAM VIBSYS” is 
presented which performs the mathematical operations of the equations 
developed in the previous chapter. The program is coded in Fortran 60 



tion 1604 computer. Although the Fortran 60 language is applicable to 
most large digital computers there are minor variations which are pecul- 
iar to the specific system in use. The Fortran 60 language does not per- 
mit automatic operations with complex numbers; therefore all operations 
involving complex numbers are accomplished by operating on the real and 
imaginary parts separately and then recombining. 

The eigenvalues and eigenvectors of the reduced system are found by 
a matrix iteration scheme which utilizes the direct and inverse power 
methods and matrix deflation. A mathematical subroutine MATSUB is used 
to carry out these operations. ^7,^ As presented, a maximum of twenty 
iterations will be performed using the direct power method and then a 
maximum of twenty using the inverse power method. If convergence is not 
reached with the inverse power method a print out to this effect will be 
executed. 

Subroutine "INVERT" is used for matrix inversion. "INVERT" uses 
the Gaussian elimination and pivotal techniques. Inversion of the 



Subroutine INVERT is a library subroutine of the Naval Postgraduate 
School and is designated locally as FI NPS INVERT. 



programming language 




specifically for the Control Data Corpora- 




matrix which contains complex elements is achieved in 
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the following manner. Let |a + be the matrix to be inverted and 
[c + jl^ the inverse to be determined. Then by definition 

[A+j6][c+jT5] = I 

The above matrix equation leads to two simultaneous equations with 



two unknowns, C and 


D. 


(a) 


AC-BD = I 


(b) 


■8C +AD - o 


Solving first for C 




C = 


[a + -ba'b] 



and then for D 

D = - a‘t3C 



Thus the complex inversion is reduced to multiplications, additions, 
and inversions, of matrices of real numbers. 

Flexibility and utility are the principal aims of the program. 

Usage requires a knowledge of the mass, stiffness, and damping matrices. 
Although various output options are available which require additional 
input data, the three above mentioned matrices are all that are necessary 
to determine the eigenvalues and eigenvectors of the reduced system and 
the natural frequencies and mode shapes of the original system, 

3,2 Program Options , In addition to finding the natural frequencies and 
mode shapes of the system, five options are available which describe the 
time behavior of the system under conditions of free and forced vibration, 
(a) Option 1, The time solution of the free vibration problem 
is obtained in general form and no additional input data is required for 
execution of this option. In section 2,2,2 the general solution of the 
homogeneous equation was developed and is repeated here for convenience. 
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(21) X (i:) [SrMi^rto) v (^XrGr^A+ /Srrt^w\ ■»-SirC)xco]e'^’"cos|5^i: 



r^i 

Kl 



^ (-o<r -/3^<5rH + H^C) X^)] 

t-il 



-(/A V 

<2 SitO/Sr'C 



The output of option 1 consists of the four coefficient matrices of the 
X(o)Co5fJr^^ XCo)Cos^A ) XC«) SINl|(Jr^ and Xto) 

terms. Therefore, the output of option 1 will consist of 4N square matrices. 

(b) Option 2. The execution of option 2 requires as additional 
input data the values of the initial displacement, and initial velocity 
vectors. The product of the coefficient matrices of option 1 and the 
initial displacement vector or initial velocity vector, as appropriate, is 
performed to obtain a coefficient column for the cosine and sine terms. 

(c) Option 4.* The general steady state solution Of the forced 

i I 

vibration problem Is provided by option 4. The general solution with a 
forcing function f(t) was developed In section 2.3.1 and found to be 



(28) 




t 

C'C) ^COS^rCV-'t) dr 

o 




f (r')e 






In this case the coefficient matrices of the convolution integral are 
evaluated. There will be 2N square matrices, N corresponding to the 
coefficient of the integral involving the cosine term and N coefficient 
matrices of the integral involving the sine term. 

(d) Option 5. The steady state solution of the special case 
of forced vibration with sinusoidal excitation is determined in option 5. 



* 

Option 3 is described in subsection (e), following. 
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In section 2.4.1 the expression for the time behavior was found to be 

(39) (ad^- b'di') coscoV - (.Bt>R+ADl) SIMCO^ 

Additional input data required for the execution of this option include 
the driving force amplitude and excitation frequency. The output consists 
of the coefficient column vectors of the cosine and sine terms. 

(e) Option 3. A plot of displacement versus time is made for 
each mass, for one of the three following cases. 

1. Free Vibration 

2. Steady state vibration with sinusoidal excitation 

3. Transient plus steady state with sinusoidal excita- 
tion 

Case 1 is obtained when the initial conditions of displacement and veloc- 
ity are given as input data and the driving force amplitudes and excita- 
tion frequency are zero. 

Case 2 is obtained when the driving force amplitudes and excitation 
frequency are given and the initial conditions of displacement and 
velocity are zero. 

Case 3 is obtained when the initial conditions, driving force amplitudes 
and excitation frequency are all given. 

Any combination of the options may be obtained with one set of input 
data. Input data format and output control is described in detail in 
Appendix B. 

3.3 Accuracy of Method . The accuracy of the solution is contingent on 
the accuracy with which the eigenvalues and eigenvectors are determined 
and the effect of computer roundoff error. Internal checks have been 
provided so that the integrity of the solution may readily be evaluated. 
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3.3.1 Internal Checks . Prior to the determination of the eigenvalues 



and eigenvectors the trace and determinate of the matrix. A, is calcu- 
lated. A check is therefore readily available since the trace of the 
matrix A is equal to the sum of the eigenvalues of the characteristic 
matrix of A and the determinant is equal to the product of the eigen- 
values. The trace and determinant of A and the sum and product of the 
eigenvalues are included as standard output in Program VI6SYS. 

Additional checks are inherent in the solution. In Option 1 the 
summation of the coefficient matrices of the X(<d coS/^r"t terms 

must sum to the identity matrix, since the evaluation of x(t) at time 
zero must equal the initial conditions of displacement. This is easily 
seen by referring to equation 21 and considering the initial values of 
velocity to be zero. Similarly, in Option 2 the sum of the coefficient 
columns of the cosine terms must equal the initial conditions of dis- 
placement. 

In the steady state solution with sinusoidal excitation, the inversion 
of the complex matrix presents one of the greatest possible sources of 
error. Therefore when Option 5 is executed, the output includes the real 
and imaginary parts of the inverse and the product of the original matrix 
and its inverse. This product should, of course, be the identity matrix. 

3.3.2 External Checks . In order to determine the reliability of 
Program VIBSYS, sample problems of two degrees of freedom for which solu- 
tions were available in the literature were programmed. Correlation of 
results was excellent. No suitable sample problems involving more than 
two degrees of freedom were found in the literature. Therefore, in order 
to extend the range of testing, solutions of undamped systems were compar- 
ed with program solutions of similar systems with negligible damping. 
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The results were as anticipated; as the damping was decreased the natural 
frequencies approached those of the undamped system* The results of a 
sample problem are demonstrated in Appendix E* 

3.4 Program Limitations . Although the digital computer extends the size 
of system for which solutions are obtainable, the ultimate size is dictat- 
ed by the storage capacity of the computer. The CDC 1604 computer has 
approximately 32,000 storage locations; however, it is easily seen that 
this may be rapidly exhausted. Program VIBSYS is limited to a system 
with 10 degrees of freedom. 

The basic principle in the development of the general equations is 
the orthogonality relations of the eigenvectors and the parameters of 
the system. Therefore each eigenvalue must have associated with it a 
unique eigenvector. This requirement is not fulfilled in the case of 
a system having two equal eigenvalues, and therefore two natural fre- 
quencies of equal magnitude. In such a case VIBSYS determines the 
frequencies and the program terminates at that point. 

If the damping is critical in any part of the system an attempt 
will be made to obtain the eigenvalues and eigenvectors. However, if 
they are found, they will not occur in the complex conjugate pairs and 
the program will be terminated. 
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CHAPTER IV 



CONCLUSIONS AND RECOMMENDATIONS 

Program VIBSYS provides an accurate method for determining the natural 
frequencies, mode shapes, and time behavior of a subcritically damped 
dynamic system. The compiling time for the program is 4 minutes and 30 
seconds, while the running time for a system with four degrees of freedom 
is approximately one minute. (Running time will depend on the number of 
iterations required to obtain the eigenvalues and eigenvectors.) There** 
fore it is more efficient, with respect to computer time, to make multiple 
runs. 

Comparison of natural frequencies of undamped and lightly damped 
systems has shown that the undamped approximation is very good for damping 
ratios of less than 0.01. For systems tested, up to four degrees of 
freedom, the natural frequencies of the damped and undamped systems did 
not differ in the first three significant figures. 

The program may be extended in a variety of possible ways. The 
general nature of Program VIBSYS requires liberal use of computer stor** 
age space. The size of the system may be extended by segmenting the 
program into separate programs for handling specific problems such as, 
free vibration, forced vibration with a general forcing function, or 
forced vibration with sinusoidal excitation. Furthermore, with minor 
modifications the present program may be terminated upon finding the 
natural frequencies and mode shapes and used as a separate program to 
obtain input data for specialized programs. 

Program VIBSYS can be augmented and made more useful by eliminating 
the restrictions which cause the program to be terminated if either (a) 
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the system has two eigenvalues of equal magnitude, or (b) one of the modes 
has aperiodic free motion, that is a purely real eigenvalue# 
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APPENDIX A 



PROGRAM STRUCTURE 

A.l General Remarks , Program VIBSYS is composed of a main body which is 
divided into five sections, and two subroutines, INVERT and MATSUB* Func- 
tion subroutines used include the sine, cosine, exponential, square root, 
and absolute value functions. These function subroutines are called from 
the Fortran library tape. Subroutine draw is also available on the For- 
tran library type and is programmed specifically for the CDC 1604 computer 
for use with the CALCOMP plotter. 

Subroutine MATSUB, which is available as a CO-OP Library mathematical 
subroutine, required minor modifications in order to retain the eigen- 
values and eigenvectors in storage for use in the main body of the program. 
The flow chart provided with MATSUB in the CO-OP Library was found to be 
inadequate for understanding. In order to perform the necessary modifica- 
tions a revised flow chart was made up and included in Appendix C. 

Main Body « The main body consists of five sections with functions 
as listed below 

(a) Input - In addition to allocation of storage spaces for 
dimensioned arrays and setting up a control for reading input data, the 
input section contains a control sequence for MATSUB and curve labeling, 

(b) Natural frequencies and mode shapes - The reduced character- 
istic matrix is constructed and the eigenvalues and eigenvectors are deter- 
mined. The natural frequencies and mode shapes are then found. 

(c) Forced Vibration - The coefficient column vectors of the 
sine and cosine terms of option 5 are calculated. 

(d) Free Vibration - The outputs of options 1, 2, and 3 are 
formulated. 
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(e) Forced Vibration with General Excitation - The coefficient 



matrices of option 4 are calculated. 



A. 3 Subroutines . 

(a) MATSUB - Evaluates the eigenvalues and eigenvectors of the 
reduced system. 

(b) INVERT - Inverts a real matrix. 
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PROGRAM NOTATION 



AB 

ABC 

ALIS 

ALRS 

AM 

AMD 

AMM 

AMS 

ARN 

ARY 

BRN 

CD 

CDC 

CPD 

DI 

DISPI 

DISPR 

DM 

DR 

EPl 

EP2 



coefficient column vector of cos/fi,.t term, output of 
option 2 

coefficient column vector of cos/3^ t term when X(0), and 

X(0) are modified to obtain total solution 

control parameter for MATSUB 

control parameter for MATSUB 

inverse of mass matrix 

product of inverse of mass matrix and damping matrix 
mass matrix Ib-sec^/in 

product of inverse of mass matrix and stiffness matrix 

-2 ©^^(V^’^MV " W^'^MW,) -4 i V^'^MW + V^'^CV - W ’’’^CW 
r r r r' | r r r r r r 

imaginary part of identity matrix 
2 (Vy''^MVj. - - 4 + ZV^'^CV^. 

coefficient column vector of sin term, output of 

option 2 

coefficient column vector of sin/^^t term when X(0) and 

X(0) are modified to obtain total solution 

coefficient column of sincot, output of option 5 

imaginary part of driving force amplitude 

imaginary part of modal matrix 

real part of modal matrix 

damping matrix lb«sec/in 

real part of driving , force amplitude 

iteration parameter for MATSUB 

iteration parameter for MATSUB 



GBI 

GBR 

GRM 

HRM 

IC 

IDET 

lEG 

IFC 

ITITLE 

IVEC 

MIT 

HITS 

MPl 

MP2 

MP3 

MP4 

MP5 

N 

RIN 

RMI 

RPI 

S 

SPEC 

STEP 

TIPI 



iteration parameter for MATSUB 
iteration parameter for MATSUB 

coefficient matrix of X(0) cos/8^C, output of option 1 

coefficient matrix of X(0) sin /Qj,t, output of option 1 

control parameter for VIBSYS 

control parameter for MATSUB 

control parameter for MATSUB 

control parameter for VIBSYS 

title for graphical output 

control parameter for MATSUB 

controls number of iterations in MATSUB, power method 
controls number of iterations in MATSUB inverse power 
method 

option 1 control 

option 2 control 

option 3 control 

option 4 control 

option 5 control 

order of system 

real part of identity matrix 

coefficient column vector of cos CJ t term, output of 
option 3 

real part of the inverse of [k- + j to 

stiffness matrix Ib/in 
spectral matrix 

step size for graphical output 

imaginary part of Inverse of K- to^ + j o C 
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UI 


- 


matrix of reduced system, imaginary part 


UR 


- 


matrix of reduced system, real part 


VALI 


- 


imaginary part of eigenvalue of reduced system 


VALR 


- 


real part of eigenvalue of reduced system C 


VCO 


• 


Initial velocity vector modified to account for initial 
velocity of steady state solution 


VCV 


- 


v/cv. 


VCW 


- 


v/cWr 


VECI 


- 


imaginary part of eigenvector of reduced system (V^) 


VECR 


- 


real part of eigenvector of reduced system (W^) 


VO 


- 


initial velocity vector 


VMV 


- 




VMW 


- 


VrV 


WCW 


- 


“rH 


WDM 


- 


imaginary part of [k - 03^ + j CO 


WMK 


- 


real part of [k - co^ + 


WMW 


- 




WO 


- 


excitation frequency 


X 


- 


time coordinate for graphs 


XCO 




initial condition of displacement modified to account for 
initial displacement of steady state solution 


xo 


- 


initial displacement vector 


xoc 


- 


coefficient matrix of X(0) cos^^-t term, output of option 


xos 


- 


coefficient matrix of X(0) sin/5j-t term, output of option 


y 


- 


ordinate of graphical output 



The above notation lists the principal array and parameter names used in 
the main body of the program. Array names not listed are used for inter' 
mediate operations. Where appropriate the names are associated with the 



symbols used in the mathematical analysis 
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APPENDIX B 



INSTRUCTIONS FOR USE OF PROGRAM 

B.l Purpose . The purpose of Program VIBSYS is to determine the natural 
frequencies, mode shapes, and time behavior of a subcritically damped, 
linear dynamic system with N degrees of freedom. (N ^ 10) The mass, 
damping and stiffness matrices of the system must be known. 

B.2 Input Data . A blank card must follow the last end card of the pro- 
gram deck. The data cards then follow in the order and format listed be 
low. The first data card uses the input format 915 and is the control 
card for the program. The nine fields are designated as follows: 



• 

11 


the order of the system 


2. MPl = 


1, Option 1 executed 
0, Omit 


3. MP2 = 


1, Option 2 executed 
0, Omit 


4. MP3 = 


1, Option 3 executed 
0, Omit 


5. MP4 = 


Option 4 executed 
0, Omit 


6. MP5 = 


1, Option 5 executed 
0, Omit 


7. IC 


1, Initial conditions given 
0, Initial conditions are zero 


8. IFC = 


1, Amplitude and frequency of excitation given 
0, Omit 


9. NPTS 


Number of points to be calculated for graphical 
output 



All nine fields must be right justified. 
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The control card is followed by the mass, damping, and stiffness 
matrices respectively. Each matrix is read in rowwise using input format 
8F10.3* with each row starting on a new card for systems with N ^ 8. For 
systems with N > 8 the elements are read in rowwise in a sequential manner, 
that is the first element of the second row should appear in the field ad- 
jacent to the last element of the first row* In this case each matrix must 
start on a new card* 

The initial displacement vector, initial velocity vector, real part 
of the driving force amplitude vector, and imaginary part of the driving 
force amplitude vector, follow, in that order, using format 8F10.3** 

Each vector must start on a new card* 

The final data card uses a 2F10.3* format. The first field contains 
the excitation frequency and the second the step size (i.e., time increment) 
for the graphical output option* 

Blank cards must be used for initial conditions of displacement and 
velocity, real and imaginary parts of forcing function amplitudes and 
final data card when these values are zero* Multiple runs may be process- 
ed by placing additional data decks behind the first* The program is 
terminated by placing a blank card behind the last data card* 

•k 

Should this format present undesirable restrictions on the size of 
the elements of the matrix the so called **E** field may be used by changing 
card number 32* The ’*E*’ field was avoided since it is more susceptible to 
error by the user* 
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B.3 Deck Assembly . The first card of the program deck is a job 
card. The statement "Use scratch tape" must be included in addition to 
the' Required job card information. 

JOB CARD 
“program VIBSYS 

END 

"subroutine invert 

END 

Tubroutine MATSUB 

END 

END 

BLANK CARD 
DATA CARDS 

BLANK CARD 

B.4 Cautions to User . The curves drawn in the graphical output option 
are straight line approximations between computed values. The step size 
must therefore be chosen appropriately in order to obtain a smooth curve. 
Since the maximum number of points permitted by subroutine DRAW is 900, it 
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may not be possible to see the transient phase completely die out. 

If difficulty is encountered in obtaining the eigenvalues and eigen- 
vectors, the control parameters of MATSUB may be altered to enhance con- 
vergence. The value of MIT and HITS, card numbers 20 and 21, control the 
maximum number of iterations to be performed in the power method and in- 
verse power method respectively. ALRS and ALIS, card numbers 78 and 79, 
represent a complex parameter that may be called the ^origin”. MATSUB 
will usually converge on the eigenvalue most distant from the ^origin**. 
ALRS = 1.0 and ALIS =0.1 in Program VIBSYS. The convergence criteria 
for the power method is set by the value of EPl, card number 85 while 
the inverse power method is determined by EP2 card number 86. EPl = 10*^ 
and EP2 = 10" in the present version. Furthermore, if lEG, card number 
76, is set equal to one, the eigenvalue iterants will be printed out and 
more suitable values of ALRS and ALIS may possibly be found by inspection 
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APPENDIX C 



FLOW CHARTS 



The symbols used in the flow charts are defined as follows: 



- arithmetic, read or print operations 




- connectors 




- decision 




call subroutine indicated 



The notation used in the flow chart for MATSUB corresponds to that used 
in the CO OP Library version, reference 7. 
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FLOW CHART M^ilH K)DY PROGRAM VIBSYS 




C-2 




C-3 






C-4 




C-5 





XCO = XO - Ri'IE 
VCO * VO - WO*CPD 



710 



ABC = GRM*VCO + XCC-^XCO 
GDC -HI2M*VCO + X03*XC0 



711 
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= ABC 
: nnn 


. € 




COIiTIMUE 




YSCALE liAX ORD/2 



STORE GliAPH TITLES 
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SUBROUTIHE MTSUB 




ITERATION Pa^ER METHOD 
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" X- 


,Y^-. 
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1 
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X-, -NORMALIZED Y; 




t>i.T FLa.1 
ICT= IQQQ 
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ISL 


= 1 • 
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VALR(N)= ALR 
VALI(H)= ALT 
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PRINT N,ALR,ALI,ICT,IJ 
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JL 
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i£.s3:...s!- soRij . 
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VALR(N) = AUR 
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VAU(N) = ALI 
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MATRIX DEFINITION 


t 402 






1 CX; =r 0<- + 


B = DEFLATED A 








IT = IH- 1 


> 




• V 
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APPENDIX E 



SAMPLE PROBLEM 

The sample problem presented is that of a four degree of freedom 
system shown in the figure below* In order to show a comparison of data 
for damped and undamped systems, the example chosen is one in which the 
natural frequencies of the undamped system may be obtained with relative 
ease* The program is not restricted to problems with the symmetry dis- 
played in the example. 




Let = M 2 - = 10 lb sec^/in 

Kj^ = K 2 = K 3 = = 1000 Ib/in 

~ 0 ^ “ “ 30 lb sec/in 

For the undamped case the frequency determinate then becomes 



2K - M ^ 

-K 

0 

0 



-K 

2K - Moi^ 
-K 
0 



0 

-K 

2K - MtO^ 

-K 



0 

0 

-K 

2K - MOO^ 



E-1 




M 



I 



o 



Expanding the determinate 

co®- 8 (^)a:»^ + + 5(1^) = 

The roots of the frequency equation are: 

CO, = 6.18034,CO^= 11.75571, CO 3 = 16.18034 , 604 = 19.02113 

The problem was programmed with a forcing function f(t) = lOOsln 
(60t) , and an Initial displacement on M^of 0.5 Inches. A step size of 
0.005 seconds was used for the graphical option. All five options were 
called for and the resulting output Is shown on the following pages. 

In addition to the output presented for C = 30 Ibsec/ln runs were made 
with C = 0.01, 0.10, and 1.0 lb sec/ln. The frequencies and damping 
ratios for all runs are tabulated below. 



C lb sec/ln 




CO, , rad /sec 




63 ^ , rad /sec 


0.01 


— 


6.18034 


0.0001 


11.75571 


0.10 


0.0003 


6.18034 


0.0006 


11.75570 


1.00 


0.0031 


6.18028 


0.0067 


11.75530 


30.00 


0.0098 


6.12699 


0.1821 


11.38430 


C lb sec/ln 


5 , 


^ rad/sec 




rad /sec 


0.01 


0.0001 


16.18034 


0.0001 


19.02113 


0.10 


0.0008 


16.18033 


0.0009 


19.02111 


1.00 


0.0081 


16.17928 


0.0095 


19.01941 


30.00 


0.2584 


15.19737 


0.3118 


17.40396 



E-2 



MATr^lX ANALYSIS OF A MJLTI*DHbRbfc OF FREEDOM VIBRATION SYSTEM WITH VISCOUS DAMPING. 
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iHb FREQUENCY OF THE FORCING FUNCTION IS 60 . 0 0 OOR AD 1 ANS PER SECOND 




o 



o 



o 

C3 

iTi 



=> r\ 

o x> 
C 3 o 
O O 

O o 

o ir\ 



rsj 



t 

Ul 

r 

<u 

7 

2 

\ 

u 

u 

T 

w 

<u 

> 

z 



q: 

s, 

o 

V 



»/> 

7 

O 

or 

u 

t- 



*2 



0 

o 

IL 



^0 

0 ^ 

to 

(kf 



VO 

7 

9 

i- 



4 

'•u 

V- 



J 

a 

y 

& 



jj 

■JO 

o 

r-o 

AJ 

T) 



(\i 



UJ 

CO 

o 

ro 

f\j 

00 



CM 

o 

♦ 

UJ 



LT\ 

r\ 

<o 

o 

o 

IS 



CM 

o 



rs. 

r\ 



CM 



LU 

M- 

O 



fN. 

IS 



I 

UJ 



t ^ W »-• O O «H 
I O O O t— I »-t r-4 



^ ^ »-f O O 

O O O O O O < 



> IS lf» 

» o O 

> •>. (/I 

> CM ^ 

> O 

Ct 

I • o 

• I t— 



J UJ LU UJ 

• M* 

) CM rsj ko 

./N m M- 
^ M’ m 
uf> m vr 
o o ^ 
rv 

ro »o m 



LO UJ 
OCM 



Ui Ui < 
OCM 



I 



I 



LLI UJ UJ LLI LL. UJ 

ro o o ro in m 

X rH OD O CO 

ro X X ro o o 

,H rv w X ro 

M- X X M- o rv 

»H ro ro »H o 

OCM M-mr>.i/>t-rc5o^oo 

— o m m ro th 



I 



CMro t 
o o < 
^ ro i 
o m 



5 tH in c> CM ' 



O fV. o 

O M" tH 

o m m 



ooocoooo 

I ♦ I ^ 

III III 111 1)1 i; 111 1 1 1 Ml f-> 

rooorowooCM ♦ 

xwtHX'nr».or«vtjj 

roxxro^CMoroo 

f-trx.fN.THM’ rx.ooTH 

M xxM•rx.®oc^Ln 
tHrorOfHOooCMo 
I O «H tH M’ CD tH 



o O rs. 

OCM CM^rommrooM- oth 



rs. CO O 

CM « m 

o a: • 

• o 

I r— 



o o o o 
till 
UJ UJ UJ UJ 
M- rs rs M 
X CM CM X 
M- in m M 

in M M- m 
in j> in if» 
o o c^ 
rs rs c> 
ro ro in 

I 



O O CD CD 
fH »H tH O 
I I I ♦ fH 
UJ UJ UJ UJ o 
CM rs. ro CD ♦ 
O O O O UJ 

X o ro o CD 
fs. CM O O O 
OM-OOM- 
tH tH rs o o 
o tH C> o CM 
th rs ro o rs 



> o o o o 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
<3 O O O O CM 
CM O O CM CD m 
tH CM CM tH O X 
CM O O CM O X 
ro M- M- ro o o 
rs rs. rxv rs o tH 

to ro tH tH ro o ro 

— • X in in X o M- 



I I 

UJ UJ 

ro ro 
o o 

rs rs 

th ro 

M rs 

c^ o 



fHfHfHtHOX OX 

OOOOOO OO 

fH I I I I 4 - t t I 

O UJ UJ UJ U.I UJ UJ UJ UJ 

^ oooooM* inx 

UJ fMOOCMoO OO 

o iHCMcmthox rsrs 

O CMOOCMOtH ox 

rf rOM-M-rOOM- Xm 

o rsrsrsrsoM no 

CMtOrOfHfHrO OO J>0 

rs. — •Xinn X oro rscM 



CM OC I 

• o 



•I UJ C 3 o 
X z> — I 
_r UI u 
T <I O 
»— > C 5 O 



'.5 xat o 
— « — o 
UJ L> Xl 

o w 

I X '/) • 

I ►- CO I 



CM Cm tH 
o o o 

I I I 

UI UI UI 
X X o 
o o o 
O O o 
in in o 
•3 c 3 O 
fs. rs o 
CM CM in 
O O fH 



fHOOfH II UJCMfH 
OO OOUiCjOO 

♦ ♦ ♦ ♦ — I I 

ujujuja._iuju'uj 
X O 
Q O O 
UJ o o 
I— n o 
■•too 
« rx o 
c_) CM in 
O O fH 



o o o o < 
o O O 
oro ro 
o ro o 



o X 

O fH 
tH <3 



o o c o 
I I ♦ ♦ 

UJ UI LL UJ 
0 X 00 

o o o o 
CD o ro o 
o in ro o 
o o o 
o rs X 
in CM th 
tH O X 



tH O II 
CD O U 
♦ ♦ ♦ 3 

UJ UJ UJ _l 

o o o < 

0 0 0 3 

o O ro 2T 
O o ro UJ 
o o o o 
o o X — • 

O O tH ‘UJ 
tH tH X 



UJ CM tH tH CM 

C!3 O O O 



JJ UJ UJ UJ UJ UJ I 
X O O X O I 
Q O O O O O I 



O fH tH O II JJ fH C 
JJ CD CD I 
Z) ^ I 
-J UJ JJ 



I tH ,H O O tH II 



O O O o 

* * * * 



UJ I 



> ro I 



in o o in ro < 

•4000001 

^ fs. c o rs X I 
CO CM in in CM th I 
O O tH fH o X < 



LL* UJ UJ 

< O X X o 

‘ 3 Q O O O o 

I 2? UJ O O O O 

I ijj r— o in in o 
> O •< o o o o 
— • CD rs rs o 

I JJ o in A) CM in 

» O tH o O tH 



■ I — n I 



/) I 



I 



I to 
•- J) 



o o o o JJ CD 

♦ * * ♦ ID — 

Ur JJ UJ UJ J UJ 

o o O o < 

o o o o > a 

o ro ro o 2 T 'JJ 

o ro ro cDjji- 
o o o o CD 

o X X o ^ 

O fH fH O JJ 

tH X X fH 



fH CM CM tH tH O O 5 
O O O C O O O I 
^ I I I I 4 ^ -f -f 

UIUJUJUJUJUJUJUJ. 

o in X o o o o I 
QOOOOOOOi 
'jjooooororo i 
I— oinmoororo i 

400000001 
— forsrscDoxxi 
cj in CM CM in otHth < 
OfHOOfHfHXX’ 



I M JJ tH CM CM T-I 

> JJ LD O O O O 

ID — I I I 

I J UJ UI LU UJ UJ 

>4 O O O O 
» 3 O CD O O O 

> 2 T JJ O tH O O 

> UJ I— o n in CD ( 

> CD -<0000 

I M ^ c rs rs c 

> UJ CD n Aj CM in 

I O CD' O 



♦ * 

UI LU 
O O 

o o 

OlO 



ro o 
ro o 
O o 



I 



II — 1/5 I 



X 



ITl 

CM 



ON 



♦ 

Ui 



UJ 


CM 


o 




w-i O 


Ok 


O O Ok 


tH *-• w cr> 


o o o o 


3 


O 


o 


O 


O O 


o 


o o o 


»4 iH 


o o o o 


-J 


4 


4 


4 


4 4 


1 


• 4 1 


• lit 


4 4 4 4 


■< 


LU 


LU 


LU 


LU LU 


LU 


LU LU LU 


LU LU LU LU 


LU LU LU UJ 


> 


lA 




O 


o Ok 


«4 


O OCM 


o CO 


<3000 


Z 


m 


O 


O 


o Ok 


LA 


fv O fk. 


fkk oro ro 


M> 00 00 M> 


LU 


<> 


fO 


o 


o ro 


w4 


CM o ro 


kO fk. fv ro 


kO kO kO kO 


C9 


ro 


ro 


o 


o ro 


V 


rv o o 


Ok kC> kO kO 


Ok o o o 




o 


o 


o 


o o 




CO o o 


v4CM kO kO 


kO «-t 


LU 




flO 


o 


O CO 


o 


o o CM 


tA lA O 00 


00 lo ro fio 




rv 




o 


O tH 


•4 


O ^ 


CM LA O 


roCM CM ro 


U. 






wi 




kO 


M* o ^ 


ro kO Ok 


*4 CM CM ^ 


O 


• 


• 


• 


• • 


• 


• • • 


• • • • 


• • • • 



u 

o 

o 

0 £ 





































CL 




CM 




kH 


<3 


O 


kH 




<3 


o o o 




tH kH tH kH 




o o o o 








































O 




o 


O 


o 


O 




kH 


kH kH <3 




kH kH kH kH 




o o o o 








































4 




4 


4 


4 


4 




1 


1 1 4 




• III 




4 4 4 4 




































C> 




LU 




LU 


LU 


LU 


LU 




LU 


LU LU LU 




LU LU LU UJ 




LU LU LU LU 


kH 
















kH 


















<3 








O 


Ok 


Ok 


O 




CM 


fs. ro o 




Ok Ok 4 O 




o o o o 


O 
















O 


















• 








o 


c> 


c> 


<3 




45 


Ok C^ C3 




45 kH kH lO 




004 4 00 


♦ 
















4 


















LU 




ro 




o 


ro 


ro 


O 




00 


4) ro o 




ro oo Ok rv 




45 ro ro 4> 


LU 
















LU 


















kH 




rv 




o 


ro 


ro 


<3 




rv 


CM o o 




kO 45 kH 45 




o o o o 


O 
















Ok 






















Ok 




o 


o 


o 


O 




o 


HT C^ o 




ro LA CM 4> 




45 00 oo 45 


















M- 


















« 




kH 




o 


oo 


oo 


<3 




kH 


kH fv o 




kHCM O O 




ro kH kH ro 


CM 
















CM 


















rv 




LA 




o 


kH 


kH 


O 




Ok 


kH Ok CD 




C^ ^ O 4 




CM 45 4> CM 


fkk 
















rv 


















o 




kH 




kH 


O 


45 


kH 




kH 


ro o 




CMCM kH Ok 




CMro ro CM 


ro 
















ro 


















CM 




• 




• 


• 


• 


• 




• 


• • • 


r— 


• • • • 




• • » • 


A 
















LA 


















OO 










1 


1 






1 


1 


4 


1 • 






rH 


































LA 






















LA 








O 
















kO 


















• 






















kH 








• 
















• 








































• 




»— 




1 
























































4> 




4 

LA 

tH 








»4 


o 


O t-i 


O 


O 


o 


c> 




kH 


o 


o 


kH 


O 


o 


O Ok 






CM 




kH 


O 


O 


kH 




<3 


00 o oo 




kH kH kH kH 


• 


o o o o 






O 


o 


O o 


kH 


o 


kH 


O 




O 


o 


o 


o 


kH 


o 


kH OCM 




O 




O 


o 


o 


O 




O 


O O O 


w 


kH kH kH kH 




o o o o 


O 




1 


4 


4 1 


1 


4 


1 


1 O 




1 


4 


4 


1 


1 


4 


1 1 


o 




4 




4 


4 


4 


4 




4 


1 1 1 


in 


III! 




4 4 4 4 


O 




LU 


LU 


LU LU 


LU 


LU lU 


LU O 




LU 


LU 


LU 


LU 


LU 


LU 


LU LU 


4 




LU 




LU 


LU 


LU 


LU 




LU 


LU LU LU 


O LU LU LU LU 




LU LU LU LU 


♦ 




t4 




kH kH 


V 


o 


kH 


O 4 




kH 


kH 


kH 


kH 


C^ 


O 


tk. 


LU 




ro 




O 


O 


Ok 


O 




O 


4” lA oo 


O LA CM Ok 4 




o o o o 


LU 




O 




kO kO 


00 


O LA 


ro LU 




>o 




kO 


yO ro 


<3 


Ok 00 


O 




o 




o 


Ok 


Ok 


<3 




O 


ro Ok o 




C>tA kH CD 


w 


004 4 00 


CM 




ro 


O 


oro 


kH 


O LA 


kH CM 




ro 




kO ro 


OO 


O 


o ro 


O 




ro 




o 


ro 


ro 


O 




o 


4) fv rv 


X ro LA oo rv 


CO 


45 ro #o 45 


O 




Ok 


o 


O Ok 


«H 


o ro 


CO O 




c> 


O 


o 


Ok LA 


O 


CM Ok 


o 








o 


ro 


ro 


O 




o 


kH Ok 00 




45 LA 45 45 


o 


O O O CD 


<> 




kO 


iH 


kHO LA 


o 


ro 


LA O' 




kO 


kH 


tH kO LA 


o 


LA 


o 




OO 




o 


<3 


o 


o 




o 


4* 4> LA 


-1 rv ro LA CM 


L5 


45 00 00 4) 


M- 




LA 


tH 


kH lA 


O 


O CM 


CM M- 




LA 


kH 


th la 


Ok 


o 


LA fkk. 


o 




ro 




o 


00 


00 


o 




o 


4- LA O 


LU 


454 CM LA 




ro «H kH ro 


Ok 


(/) o 


M> 


kO O 


o 


O 


■«r 


kO Ok 


to 


C> ^ 




c^ 


kH 


o 


00 00 


o 


LU 


kH 




<3 


kH 


kH 


o 




o 


OLA C^ 


> r>- CM 4) LA 


X 


CM 4) 4) CM 


CM 




Ok 


kH 


kH Ok »H 


O kO 


tH cm 




Ok tH 


kH 


Ok CM 


<3 


fkk. kH 


kT 


OC 






kH 


45 


4) 


kH 




o 


ro fv CM 




LAfv CM 4> 




CMro »0 CM 


fv 


































CM 


<< 


• 


»— 


• 


• 


• 


• 


»— 


• 


• • • 


Li. 


• • • • 


t4. 


• • • • 


LA 


(X 












t 


1 in 


q: 


r 


1 


1 


1 


> 






• 






cr 






1 


1 


X 




1 


O 


• 1 1 


O 




• 


o 














• 


o 
















» 


in 




HC 










-4 














1 


►— 














• 


4— 


















LU 




Q. 










X 






X 




X 






CJ 
















C_> 










































»>• 






LU 
















LU 


















u 














o 






X 




X 






3- 
















» 


















z 




< 










•< 






4- 




ft— 






Z 
















Z 
















M 


LU 




UJ 










X 






<4 




4 




II 


LU CM 


kH 


kH CM 


O 


kH 


kH 


O II 


LU CM 


kH 


kH 


CM 


o 


tH 


kH O (/) 


3 


kH 


X 


O 


kH 


kH 


o 




o 


CD <3 Ok 


X 


kH tH kH kH 


X 


CD <3 O <3 


LU CO 


o 


O 


O O 


o 


o 


O 


O LU 


to 


<3 


O 


O 


O 


o 


O 


O O 


LU 


O 


O 




o 


O 


O 


o 




kH 


O kH O 




kH kH kH kH 




o o o o 


3 




1 


1 


1 1 


4- 


4 


4 


4 3 


•»> 


1 


1 


1 


1 


4 


4 


4 4 


Z> 


LU 


4 




4 


4 


4 


4 


•k 


1 


4 1 • 


I— 


• 1 1 ft 


ft— 


♦ ♦44 




LU 


LU 


LU 


LU LU 


LU 


LU 


LU 


LU J LU LU 


LU 


LU 


UJ 


LU 


LU 


LU LU ^ 


q: 


LU 


X 


LU 


LU 


LU 


LU 


X 


LU 


LU UJ LU 


Z 


UJ (U UJ UJ 


z 


lU UJ LU UJ 






00 


O 


o oo 


O' 


O 


O 


C^ < 




CO 


O 


O 


00 


Ok 


O 


O Ok 


< 


Lk 


■4“ 


•— » 


Ok 


O 


O 


Ok 




O 


o 4- rk. 


UJ 


4>LA Ok fs. 


LU 


o o o o 


> 


Q Ok 


O 


O Ok 


Ok 


o 


o 


O > 


Q Ok 


o 


o 


Ok 


Ok 


o 


O O 


> 






X 


Ok 


O 


o 


o 


X 


ro 


O Ok 00 




CM Ok 45 fv 


mm 


4> 00 eo 4> 


Z 


LU 


o 


o 


O O 


ro 


o 


o 


ro Z 


LU 


o 


o 


o 


o 


ro 


o 


o ro 


z 


_jr Ok 


♦— 


ro 


o 


<3 


ro 


r~ 


00 


o Ok ro 


o 


oro ro 4> 


o 


4>4> kO 4> 


LU 


1— LA 


o 


O LA 


ro 


o 


o 


ro LU 


I— LA 


o 


o 


LA 


ro 


o 


ofO 


LU 


-< 


Ok 


-< 


ro 


o 


O 


ro 


•< 


LA 


OCM Ok 




4 « 4> C^ 




Ok 03 O O 


C3 


■< O 


O 


O O 


o 


o 


<3 


O IS 




<3 


o 


o 


O 


o 


o 


o o CO 


X 


o 


X 


o 


o 


o 


o 


X 


LA 


O *4 LA 


ll. 


Airv ro kH 


14. 


kH4) k© kH 


»— • 




fk- 


O 


o rv 


00 


o 


O 


CO — 


>-4 


rv 


o 


o 


fv 


oo 


o 


O 00 




3 


CM 




oo 


CD 


<3 


oo 




Ok 


OLA r«k 


u. 


45 45 kH lA 


Lh. 


OO ro ro OO 


LU 


O CM 


LA 


LA CM 


kH 


o 


o 


kH LU 


O CM LA LA 


CM 


kH 


o 


O kH 


LU 


r— 


kH 


-1 




o 


O 


kH 


-J 


kH 


O 00 oo 


LU 


ro Ok CM 


LU 


roCM CM ro 




O 


Ok 


kH 


kH O' 


k© kH 


kH 


<J 


O O 


kH 


kH 


Ok 




kH 


kH 




■< 


kC 


HK 


k© 


kH 


kH 


4> 




CM 


O fv kH 


o 


4IA CM to 


o 


vH CM CM kH 


X 


in 






























U. 


Z 


• 


Q 


• 


• 


• 


• 


Q 


• 


• • • 


u 


• • • • 


o 


• • • • 


k— 


tn 


1 


1 


1 1 








< 


CO 


1 


1 


t 


1 








O 






O 










O 


1 






1 1 1 








< 














Z 


<C 


















LU 




X 










X 






LU 




LU 




CM 


































X 


I 




















I 




I 




















u 


















3 


>— 




















1— 




1— 










* 



C 3 O C 3 O 
t I < I 

UJ LU LU UJ 



00 o 

CC (T) *0 
>c «o rv 
ro -o ir\ 
t>o r> 
o o OJ 
fo rvj 



CD CD CD CD 
I I i I 
•JJ Li.' UJ UJ 
<> ifN m o 

r> 't r> 

O tH y-i O' 
>C O' O' -o 

■o o 

X 00 CC 00 

rsj o o c\j 

,H rv. 00 T-» 



0 0^0 
tH »-• t-« W 
(lit 

U.< UJ UJ U- 
O' 00 m '»■ 
»-l Xt W O' 
fO 00 CD W 

Co ■'T rv Lf> 

O LTN X vO 
-X) O O' 
CO ^ o 
o ro >«■ 



CD CD CD CD 
o CD CD CD 

♦ ♦ ♦ ♦ 



'<y X crj 
rn O O ro 

CD o O CD 

OO >0 'O CD 
^ fO ro T-I 
«0 00 00 O 

r) 00 00 ro 



UJ UJ UJ UJ 



00 ’ 



■ 00 



OC fO ro CO 
ov O' O' 

O fO fO «o 
'O OVI 00 'O 
00 ro ro (NJ 
»H O' O' y~i 
rO T-l ^ to 

I I 



UJ UJ UJ UJ 
<0 tH tH <0 

^ "O ro ^ 
^ in ir> 'T 
rv rv ■or 

fH m m f-« 

CO CD CD CO 
^ CD CD -«• 

'O ^ «o 



UJ UJ UJ UJ 

fo -«■ rv m 

^ T-l ^ fO 

«o 00 rv CO 

^ O O' CD 
rv ^ (\j CO 

ro ro 00 ^ 

^ ro OO O' 
w w ^ 



o o o o 
♦ ♦ ♦ ♦ 

UJ UJ U’ UJ 

O CD CD CD 

^ OO CO ^ 
ro vO 'O fO 

CD CD CD €3 
OO 'O 'O CO 
wr«o to t -4 
O CM OO >0 
to OO 00 K) 



o o o o 
till 

UJ Ul U> UJ 
CD CM CM OO 
vO CD CD «0 
^ <0 <0 CD 

O' in m O' 

^ -or 

CD CM CM O 

-<r fo 
OO ro 



D CM 'M X 

O ^ 'C 

O' O 'O 

rs. fN. *0 

X X X 'O 

m m *o 



Ui UJ Uj Ul 

>0 n in in 

or X 

tH o o ^ 
O' -C X O' 
t-l X X tH 
C 'C 'O X 
CD 'O ^ 

CM r> CM 



O CD CD ^ 
f-l W t-* 
I I I I 
UJ Ul UJ .u 

r> 'O ^ n 

ro W 'T 
rs. o- O 
X 'O fv 
m m 'O ® 
o o* ro o 
O' th ro ^ 
n CM CM r> 



♦ 

UJ Ul UJ Ul 
O O CD CD 

X 'C 'O x> 

>o <o <o <o 

O O' O' o 

<o ^ w <o 

'O X X fO 
CM *0 "O CM 
CM tH (M 



0 o o o 

1 I I I 

UJ Ul UJ Ul 
M- ® ® ^ 
"O 'T ^ 

O ® X O' 
•o CM CM "O 
00 ^ CM 
•O O O' X 
O' «H W O 



CD CD CD 3 



m o om 
rv «o X 0 " 
n in m XI 



till 

Ul Ul Ul Ul 
O' X CM rv 
X ^ X 

X CM rx rv 
X ro X O' 
rv X CM CM 
O tv CO 00 

Xi m *o 00 

— ' O CM ■»H ■»-H 



O CD O CD 
♦ ♦ ♦ ♦ 
UJ UJ UJ Ul 
CD CD CD CD 

X X X X 
X X X X 
o ^ O' o 
X tH tH X 
ro X X "O 
(MIO *0 OO 
(M 00 



till 

LU UJ Ul UJ 
CM 00 OO CM 
O X X D 
X X X X 
in 00 CM in 



CM « 



■ -O' 'T 

> 00 



O' X X -o- 
tH X X w 



U) 



^ o o O CD 



_ .u II Ul UJ 

tn CD CD CD CD 



_l *0 O' o *o 

jj rv *o 

~> o 0 X 0 
-o X .n *o 



X o o o o 
till 
U U UJ UJ 

X X n X 

-- X X ' 3 ’ 
^ r~l O O W 

O X X o 
oi *-« X X 
X O X X 

K 3 ro *0 CD 

ro *o fo CM 



^ Ul U U III 

CJ T-l CO X CM 

CM rs. n 
X CM X ^ CM 
»H CD cr 
_J CD X X X 

UJ O' O' o 

D* T-l X X 
t-» X CM X 

O I t 



J UJ UJ Ul 

> O CD O 
1 X X X 

> X X X 
: O' o CD 

> ^ X 

> X X ro 
J ro ro CM 
J t-l *-( 00 

I I 



3 CD 3 O 
:? • I I I 

-- JJ UI Ul U’ 

cn 'T X X ^ 

•O ■3 O' r^ 
X O X X O' 
fO CM CM ro 
_l CM O O 00 

UJ 'O O O 'O 

;> o ^ ^ O 



tH O O O O 
T-l I I I I 
LU Ul UJ Ul 

O X X w 

fO tH ^ ro 

^ X O' O X 

— rv X X rv 

CO X X X X 

CD fV fv o 

X o ^ ^ o 

o^ooro-o- 

Li. ... . 

O I I 



</) I I I I 

O Ul Ui UJ UJ 

O ^ tH X ^ 

■O' tH -O' tH 
X CM tr CM CM 
CM O' ro CD 
_ixrvxtr 
JJ 3 ro ro 

> X X X fO 

tH tH CM tH 
Li. ... . 

O I I 



X o o o o 
tH 

Ul UJ Ul Ul 
0 0 3 0 
W. X X X X 

CO X X X X 

O O O' O' CD 

CJ X tH tH X 

ro X X CO 
X CM (0 (o OO 
m tH tH 00 
L*. - . • • 

O I I 



tH CM CM tH 
W o O O O 

^ I I I I 

•— UJ UJ UJ Ul 

to CM OO 00 (M 

o X X o 
X X X X X 

X 00 ro X 

_j o- tr XT t*^ 

LU CM o O CM 

>• -O' X X tf 

tH X X tH 

u. . . . . 

O ( I 



I- I lit 

JJ UJ Ul u 
i, ^ o = o 

_ O' fVJ tVr 

LJ ro X X ro 

« r>. X X Ov 

.i. X If' 

L^ fOfO O- 

u CM X X rvj 

o CM -o r-^ O'] 

tJ • • • • 



O CD o o 
r— I I I I 
U LU UJ UJ 
JJ o X X CD 

^ tJ t? tJ tJ 

U O' tH tH o 

« X O' O' X 

J. X tH .rH X 

X X X X 



' CM I 



I CM CM ' 



> 00 



5 0 0 3 0 



?* U U I 
JJ *0 tH • 
^ X (\ f 

u o ro 



r>. ro 

tH 00 
D X O' 
^ X X 
■I O to 
^ X ^ 



5: o o 
o o 
*—•♦■» 
^ UJ UJ 
UJ o o 
— -O' X 

O ro X 
« 3 o 
-.J. X X 
U tHrO 
-Jj X 00 

O ro 00 

c_> . . 



UJ UJ 

3 3 

X 

X ro 

3 3 

X X 

ro tH 

00 X 

00 ro 



X tH tH tH tH 
3 3 0 3 



■ CM 



^ UJ UJ UJ U 
X* "O t 
« X r 
o rv I 
— X r 



ro ro X 

O O' Ov 

ro ro X 

CM CM X 

fO ro 00 

UJ tH O' O' -r-l 

O ro tH tH ro 



U X 
U CM 



£ tH tH tH tH 
0 0 0 3 
»— I I I I 

? UJ UJ Ul Ul 

JJ n tH o X 
— -O' ro ro 
o-o-xx-o- 

— . t» fs^ rv tr 
U tH X X tH 
Li. X O O X 
LiJ -O’ 3 3 -O' 
O X -O' -O' X 



s: 00 tH 00 ro 

tH tH tH tH 

*— I ( • ( 

^ UJ U U UJ 

JJ rx. O' ro 

o ^ X ^ 

O O CM X X 

— XCMXM- 
Li. o X 0 '' 

Li. X CD O ro 

'JJ X X X -O' 

o ro tH O' X 

o • • • • 



< 

51 3 O 3 O 
o o o o 
^ * 
Z UJ U' UJ UJ 

Ui CD O O O 

-Htroox-o- 

O ro X X ro 

« O CD O O 

Li. X X X X 

U tH ro ro tH 

•JJ X CM CM X 

O ro CM CM ro 

O • • • • 

I I 



3 CD O CD 
r— I I ( I 
? UJ UJ UJ Ul 

UJ X CM 00 X 

ft ^ o CD X 

C.> X X X X 

— O' X X O 
U ^ ^ M- M' 

u O CM 00 O 

Ul ro -O' -O' ro 
O CM tH tH 00 

u • • • • 



UJ 

I 



I i 



THt COfcrnCItM MATRIX OF X SIN( 15 . 697 T) 



nH n-l n-t 


n-f 


ft nH nH 


nH 


o o o o 


CM nH nH CM 


nH nH 


nH nH 


O O O 


o 


nH nH nH 


nH 


3 3 0 0 


O O O O 


o o 


O O 


1 1 1 


1 


t 1 1 


1 


♦ ♦ ♦ ♦ 


» 1 t 1 


t 1 


1 1 


LU UJ UJ 


LU 


LLi LU LU 


LU 


LU LU LU LU 


LU LU LU UJ 


LU UJ 


UJ LU 


Ox XV XV 


rv 


vO nr o 


to 


3 0 3 0 


CM nr nr CM 


>o nr 


nr -o 


o o o 


o 


O C5 nr 


XV 


c X X X 


to nr nr to 


X to 


fO X 


O CM CM 


o 


nH rv nH 


X 


o o o o 


00 XV XV CM 


X 


rv X 


X to to 


X 


nr <o 00 


C3 


0 0 3 0 


X XV XV XV 


O XV 


XV o 


n-t nr nr 


nH 


nr -JV to 


XV 


nH 'O nO nH 


3 'O >o o 


ro o 


o to 


XV o t> 


XV 


>0 ft CVJ 


X 


X to ro X 


X CM 00 X 


nH XV 


XV nH 


o XV XV 


O 


X X CM 


3 


to CM 00 fO 


XV CM CM XV 


nH O 


O nH 


O XV XV 


o 


tH to o 


Ox 


nH CM 00 nH 


fxx nH nH pv 


nr 'O 


O nr 


• • • 


• 


. . • 


« 


• • • . 


• • • • 


• • 


• • 



0 o o o 

1 t t I 

OJ lit Ol J.T 

tr» rv rvj ir> 
O O O o 
rsi CO CO CM 
ro ^ V ro 
tj- r. rv 
O XV c> 
XV ^ ^ XV 
XV fo irv 



UtJ JJ lU JJ 
m O' fo C 3 
sr *0 ITv 
'v. vO c^ »H 
fvj 

^ CO c^ ro 
x> ^ (M rvj 
?si XV ^ rvj 

x» C\J »-» o 



o o o o 

+ ♦ ♦ ♦ 

LU UJ LU LU 
C3 O O O 
CO CO 

>0 ro ro sO 

C3 c: O C3 
O CO CO >o 
ro »H rH r«o 
Aj >o >o rvi 
cvj ro ro 00 



3 0 0 0 
till 
LU OJ LU LU 

Ow fs. ^ 

^ -O O ^ 
XV O O IfV 
XV >0 >0 XV 
>0 ^ ^ 'P 
00 X) CO oo 
rvi o o <\» 



I ♦ ♦ I 

UJ LU Ol LU 
^ CVJ CJ ^ 
rO O O "O 

rs. XV XV rv. 



XV < 



3 XV 



o rv o 

XV O' O' XV 
>0 C 3 C 3 O 
o ^ ^ >o 

I I 



0 0 0 3 
I I r ( 
lU LU LU XJ 

XV rvj c\j ,r» 

O O C3 o 
CM CO CO rvi 
K) M- fO 
M- ov rs. <«■ 
O' X XV O' 
XV ^ ^ XV 

X to ro XV 



.O I I I I 

O ’ii LU UJ LU 

U » o 
^ XV *0 o 
x: O' rvj <0 rv 
■V CM rn 'C 

_J O T-^ lO XV 

JJ <5 O O’ tH 

2* O CM XV I> 

^ >£> rv fo 



I 



I 



ai LU LU LU 
o o o o 

— ' rt^ 'T nr X) 

CO O »0 to vO 

c: o o o o 

U >0 3 ) oo >0 

rO »-l rH to 
X CM >o o rvj 

(AJ to to 00 

U. . . . . 

O « I 



— LU LU LU LU 
^ rv fs. 5 T 
nr >0 >0 M" 



X XV < 



> XV 



X 'O 'O X 
vO nr nr 'O 
JJ OVI ® oo CM 
CM O O CM 



»-t I ♦ ♦ I 

UJ UJ UJ LU 
nr CM 00 nr 

W -0 >o to 
zrn^xixvO' 
— . XV c O XV 
to o O' rs. >o 
XV O' O" XV 
X nC O O O 
>0 n-t r'l O 



O O O CJ 
till 
LU UJ a UJ 
Ov XV XV Ov 
O O C O 
O CM CM o 
X ro ro .-X) 
ft nr nr ft 
X O O XV 
O XV XV o 
O XV XV o 



z: oo ft nH CM 



— I I I I 

^ LU LLi 'JJ LL- 

X' <o r> XV o 

« o nr o 

o o O f" ft 

~ *0 nr VO nr 

L 4 . n-t "O rv 

u. OC <» XV >o 

JJ 3 O CM CO 

O n-t fl XV *-» 

(_)♦*•• 
I I 



>o 
u >o 
^ o 

U- n-l 
U. X 
UJ "O 
O n-l 



LU JJ LU 

0 o o 
X X 'O 
>0 O nO 

o o o 

>0 nO nH 

*0 ro X 
CM 00 to 
CM CM ft 

1 I 



I 



I 



Z LU U-I LL LU 

JJ CM nr CM 

— to nr nr ro 

O CM XV XV 00 

—• XV XV XV XV 

Lt. O >0 >0 o 

Lt. X CM CM X 

JJ XV CM CM XV 

O O' nH nH Ox 

o • • • • 

I I 



31 nH ft nH ft 

O O O O 
t— I I I I 
Z U-’ UJ LU UJ 

JJ >0 ^ nr >o 

— X to to X 

o X rn* X 

— O XV X o 

Lc to O >0 to 

U. r-t XV iXV nH 
lU n-t >0 X n-l 

o nr X >0 nr 

u • • • • 

I I 



IHt COEfriCIENT COLUMN FOR bXP( -.573r)C05C 

.6910E-01 .lllSE^OO .1118E*00 .6910E-01 



o# 

o 

OJ 

fO 

fO 

'tr 

vO 



CJ 

o 

Oj 

o 



o 

I 

UJ 



o 



o 

o 

CO 



o 

I 

LU 



o 

'O 



o 

IIJ 



o 

CL 

z 

UJ 

LU 

00 



CO 

X* 

X 



o 

UJ 



o 

_I 

CL 

z 

UJ 

UJ 

m 

CO 



o 

UJ 



UJ 

UJ 

CD 



o 

t 

UJ 



C 3 

o o 

♦ • 

^ u» --- UJ 

k- CO *— fO 

tH »H O 

t-i rv o 

lO tH lO c\j 



O tH 

C 3 O 

♦ » 

U) UJ 

*- ao ►— 

rv. T-» fv O 

O o 

vO 'O C\J 



ro 

CVJ 



C5 




I— DO 
C3 CSJ 

ro ro 
t\j ro 



<3 



'/> I tis 



CO I CO t 









r 



cn z 

o — 

O CO 



to 

o 

o 



to 



to 

o 

o 



to 



FO ro 

-H O »-f 

T O O O O O 

• I .4 ■ I 

'OJ t>J li (\J UJ 

■r-0 » oD I ro 

^ t-» O 

O -rH O 

w >.■ «H CV 

OL • CL • CL • 

X XI XI 

Ul UJ UJ 



JX fx 

C\J O CNJ »H CVJ 

O O O O ^ 

• 4 • I 

ro UJ ro UJ tfV 

I DO I rx I 

tH o 

r-i fX 

w W w CM ^ 

CL • CL • CL 

X « XI X 

Ul LU UJ 



fx 

C3 CVJ ^ 

O ^ O 

♦ • I 

Ul (TV UJ 

GO I 00 

*H CVI 

»H ro 

w fO 
• CL • 

X 

UJ 



GC 

o 



ac oc 

o o 



n: ac 

o o 



cc 

o 



z z z 

X 2 : i 

3 3 3 

_J 3 _l 

O CM O O O ^ 

U o CJ O (J C3 

I + I 

I— UJ ♦- LU I— )-L 

Z ro Z O z tH 

UJ ro UJ o UJ ^ 

— • ^ — OO »— tM 

O'© O tH O ro 




O o O tH O 

u O o C 3 CJ 



4 I 

I- UJ 4 - 1 X 1 4 - 

Z O Z vO 2 

UJ o UJ (VI UJ 

OO ^ in 

tj tH CJ ^ CJ 



tH O tH 

O CJ O 

I I 

UJ 4 - UJ 

O 2 IX 

•rt Ul m 

O — • O 

O CJ CVJ 



u. 

u. 

UJ 

o 

tj 



UJ UJ 

o o 

O CJ 



u u 

Li. U 

UJ LU 

o o 

O CJ 



I u « 



o 

CJ 



o 

LU 



CJ 

lU 

X 

UJ 



UJ 

LU 

33 

CO 

< 

I 

CVI 



o 

z 

2 

o 



CL 

o 

z 

o 



cr 

33 



o •- 

X <l 



X to 

o 



CO 

CO 
3 
• -J 

3 CL 



X. 

^ CC 



X XJ 
O I- 
X t* 
I— 
X CO 

o 

3 X 



LU 

to 

CO 
3 
- _l 

3 (L 



O 4- 
X <t 

t— 
X CO 

o 

3 >- 
O 
<f 
UJ 

CO 

CO 

3 



> LU 



3 C <* 
•— CK 



2 
• Ul 
;0 

O CO 
-I 2 

— CC 



C=J 

UJ 



I 

CL 

<* 

X 

CD 



I 

3 . 

< 

X 

-D 



bHAPH T1TL6U . . MIKLOS. T.J. JOB BOX M 

transient plus steady state . . HAS BEEN PLOTTED 



to 



«-l (V OJ CVI •r-l 
»H r —4 T-i 

I I I I 

UJ LU LD LU UJ 
Q ^ O 00 
o rv o lo fo 
s: o rv r>v to 
o <0 <c <0 
ac T-< c\j -o o 
O lf> 10 O CO 
U. C\J lA ^ O 
fO >0 C> 



-r-f (V (Nt C\J CVI 

o o o o 

till 

LU UJ UJ UJ UJ 
O O O O C 3 

o o ru oj o 
s: to CO CO to 

rv «o ^o pv 

tr lA to to irv 
o PO to ^ 

Uw AJ <0 <0 C\J 

CVJ 10 to CVJ 



cr 

UJ 



QC 

UJ 



rsj t -4 AJ ^ 




UJ UJ a' UJ UJ 
a o Aj Lr -^ 
O »H IT tH O' 

r to cvj o 
Cvj ^ Pv lA 
QC a* (A X o 
0 X 000 

U_ to X V <o 

•»r vo to vr 

s: , . . . 

q: I • 

jj 



CM CVS AJ CM CM 
0000 
I I I I 
UJ UJ UJ UJ UJ 
Q CM M* M- CM 
O X "O »0 CO 
i: pv o o p<v 
xo 'O to >0 

ac vO CM CM o 
O CM to »0 CM 

T-t O O tH 

to iH to 

s: • • • • 

cr t t 

UJ 



to 



to CM CM CM CM 



till 

UJ UJ UJ UJ UJ 

Q to ^ fv lA 

O M* -fH to 

X -o CM rv CO 

M o o o 
□c rv ^ CM ® 
otoiocM-xr 
Uw M” to CM o 

lA »H ^ 

X: . . . . 

QC I 

UJ 



I I I I 

Ul UJ UJ UJ UJ 
a X CM CM X 
O X o o X 
£ X X X X 
O lA lA O 
cr M* M* ^ M” 

O O CM CM o 

U. to ^ to 

CM »-• »-l CM 

i: . . . . 

cr I I 

UJ 



M- CM CM CM CM 
tH tH tH tH 
lilt 

UJ UJ UJ UJ UJ 

Q X M* O to 

o O o M* lA 

i: tH p<v tH X 

V X CM o 

cr M” lA to lA 

O X tH CM X 

X X CM O 
tH to O fv 

X • • • • 

cr I I i 

UJ 



^ 10 CM CM to 
0000 

I I I I 

UJ UJ UJ Ul UJ 

Q CM ^ M- CM 

o to ^ to 

X CM lA lA CM 

lA tA (A lA 
cr o X X o 

O X CM CM X 

u. lA CM CM lA 

t'v r -4 tH P«* 

I: . . . . 

QC I I 

UJ 



in 

o 

C_J 

CM CM CM CM 

JJ tH tH tH tH 

X I t i • 
»— UJ UJ UJ UJ 
0 0^0 
u. X tH ft to 
O to X O f''- 
X X tH X 

_| to lA CM X 

tH CM O O 

cr o X o M* 
O CM CM tH o 

UJ • • • • 



2 

to 

CM CM CM CM 
UJ o o o o 
X I I I I 
»- UJ UJ UJ UJ 
O O C 3 o 
U. CM «-• tH CM 
O X CM CM X 
X M- M” X 
_j fO O OfO 
<r 10 rv rv ,r> 
cr X X X X 

C 3 fO lA lA >0 

UJ • ■ - • 

2 



to 

o 

u 

tH tH tH CM 

JJ tH tH tH 1 — < 
X I I t I 
I— (JJ UJ UJ UJ 
(O X lA 
U. fO tH 
O A- ts. C 3 
r-t X NO IV. 

_J A A X X 

« X O CM O 

r O »H CM -<T 

C 3 A CM AJ NO 

UJ • • * - 

~ I 
2 



2 

CO 

CM CM CM CM 
:JJ 0000 
X I • I I 
UJ UJ UJ UJ 
■sr X X 

u. *0 tr ^ to 
o o X X o 
NO CM CM fO 
_! CM ^ ^ CM 
•< NO O O to 
cr o tH tH c> 

tH .tH tH tH 

UJ • • • • 

— I I 
2 



to 

O 

C _3 

to to CM (M 

JJ tH tH tH tH 
X I I I I 
»— JJ UJ UJ UJ 
o X CM rv 
Uw X ^ X tH 
O X CM rv PS. 

X to X o 
_» rv X CM CM 
< O Ps (M CM 
X A A to CM 
CJ 3 O CM tH tH 



2 

to 

CM to to CM 
UJ O O O O 
X I I I I 
I— Ul JJ UJ UJ 

CM CM CM CM 

Li. O X X O 

O X X X X 

A CM CM A 

^ ^ ^ ^ 

<t CM C 3 O CM 

X M- X X NT 

C 3 tH X X tH 



»— I I 

2 



C/) 

O 

CJ 



CM CM tH CM 



X I I I I 
t- UJ UJ UJ UJ 
A O to O 

U. ^ to A NT 
o rs X o tJ 
to to rv Ai 
_t Pv X O to 
< A ^ CM CM 
X CM A tH CM 
C3 A CM tH O 
tu • • • • 



2 

A 

CM CM CM CM 

>JJ O o o o 

X I I I I 

»— UJ Ul UJ UJ 

'^'PspvM' 
u. tr X X '»r 

o A O o A 

A X X A 
_» X ^ ^ X 
< CM X X CM 

X AJ C> O (M 

O tH tH tH tH 

UJ • • • • 

»— I I 

2 



CT> 

I 

uJ 



2 2 

o o 



2 2 

O O 



2 2 

O O 



2 2 

O O 



I- CM cmVm cm 

^ tH tH tH tH 
_1 I I I • 

O UJ UJ UJ UJ 
> A CM O 
2 O A tH O 
oto A X Ps 
O X A X X 
ps fO A CM 
JJ X ^ AJ A 
T ps CM X A 
I- A ps CM X 



»- CM CM CM CM 
Z 5 o O C 3 o 
J ! * I I 

O (JJ UJ UJ (JJ 

> C 3 o O O 
2 AJ tH tH CM 
O X AJ CM X 
O X '<P X 

to C> O to 

JJ NO ps ps fO 
X X X X X 
I— to A A to 



JJ 

X 




1(11 
UJ UJ (JJ UJ 
tH -O X CM 
A) Ps tH A 
CM A '•S' AJ 
tH C 3 X M” 
O X A A 
O tH O' O 
r-l X »H X 
tH X CM X 



I— CM CM CM CM 
X) o o o o 
-I I i • I 



O UJ UJ UJ UJ 
^ X x^ 
2 to -(r M- NO 

O O X X ^ 

u to AJ CM to 

CM M* M* CM 
JJ NO c^ O to 

X O tH tH o 



l~ CM CM to CM 

D tH tH tH tH 

-I I • I I 

O UJ UJ UJ UJ 

> ^ tH X 

2 ^ tH ^ tH 

O CM M' CM CM 

O CM O to cr 

X Ps X M' 
■JJ C3 to ps (O 

X X X A to 

»— tH tH CM tH 



H- CM to PO CM 

O O O C 3 O 
_J I I I I 
O UJ UJ UJ UJ 
2 CM CM CM CM 
2 o X X O 
O X X X X 
U A CM CM A 
M- ^ ^ ^ 
UJ CM O C CM 
X ^ X X ^ 
»— tH X X tH 



»- tH CM CM CM 
~A tH tH tH tH 
—I I I I I 
O UJ UJ UJ UJ 
>. X o 
2 tH A to O 

o o CM X rs 

C -3 NT CM ro X 

X tH X A 
JJ X C 3 tr tH 

X O CM A X 

I— tH X CM to 



»- CM CM CM CM 
00000 
-I I • I • 
O UJ UJ LU Ul 
> ♦ rs Ps. ♦ 
2 ^ X X ^ 
o A O O A 
C _5 A X X A 
X ^ ^ X 
'JJ CM X X CM 
X CM O C> CM 



Li. I I I Lis 

O O 



U. I I U. I I 

o o 



X 



X 



X 



X 



X 



X 



X 



X 



cr 






or 



X 



X 



X 



X 



X 



< CM CM CM CM 
t tH tH tH tH 
till 

»- LU UJ UJ UJ 
2 X A C> ps 
JJ CM O X ps 

— Oto to X 

tJ M* X X l> 
CM ps to tH 
Li. X X tH A 
.i. to ps O CM 
LU ^ A CM to 
O • • • • 

(Jill 



■a CM CM (M CM 
X O O O O 
I I I I 
»— UJ U U-' LU 
2 0 0 0 0 
JJ O' CM AJ C> 
—* to X X NO 

u rs X X ps 

—• A to to A 

U. to NO 
u. CM X X CM 
UJ CM to to CM 

o - • • • 

o 



<C O tH tH tH 
X tH tH tH tH 
1)11 
l~ uu LU LU UJ 

2 NO tH to l> 

UJ X Al to tH 
« O CM PS fO 
CJ ^ tH tH CM 
— • fs o A O 
Li. O X X 
Li. to tH O to 
UJ tH tH A ^ 
O • » ' • 



'I CM CM CM CM 
X o o o o 
I I I I 
I— UJ (Jj LU LU 
2 NO ^ CM 
LU X to to X 
ps O O ps 

LJ X to to X 
— X CM CM X 
u. CM to to CM 

U. tH O O tH 

X to tH tH to 

o • • • • 



< to CM to M- 

X tH tH tH tH 



lilt 

I— LU LU LU UJ 

2 ps O to 

UJ O M' X ^ 

— » O CM X X 

O A CM X 

— • O X ps ps 

U. X o o ro 

X X X A ^ 

LU to tH O A 

O • • • • 

O I I 



<* CM CM CM CM 

X «3 C 3 C 3 O 

I I I I 
*— LU X X X 

2 X CM CM X 

X X O O X 

X X X X 
CJ O A A o 



X O CM CM O 

X to Mr ^ to 

X CM tH tH CM 

O • • • • 

U I I 



< tH tH CM CM 
X tH tH tH tH 



to CM CM to 
X O C 3 o CD 



I I I I 

I— X X X X 
2 X X A X 
X O tH M* O 

O C> ps tH 
U *0 HT to ^ 
tH X ps -M- 
X X X A X 
X O CDCM X 
X tH tH A tH 

o • . • • 



I I I I 

X X X X 
2 CM 'T M- CM 

X to ^ ^ «o 
— AJ A A (M 
C_> A A A A 
— . o X X O 
I*. X CM CM X 
X A CM CM A 

X ps tH tH ps 

O • • • • 



X 

X 



X 

X 



X 

X 



POKCfcD VIBRATION OPTION NO. 4 HAS QEtN EXECJIED 



iHt HEAL. PART Of- THE INVERSE 

-.2900E-04 .5028E-06 .8121E-07 -.56276-08 

.5028E-06 -.2892E-04 .4971E-Q6 .8121E-07 

.8121E-07 .49716-06 -.2892E-Q4 .50286-06 

-.56276-08 ,81216-07 .50286-06 -.2900E-04 



NO 



Ui 

CM 

O 

iC\ 



NO 



CM 

O 

CM 



I I 



CO (T» ir\ 
o o o c. 

til' 

LiJ LU UJ UJ 
■j* CM o rv 
CM CD rv ^ 

o rN* o 

CM O r-» PO 



^ CM CM 

I I T— ( CD 
I I t ^ 
^ O! OJ UJ 
ir> r'' r-i o 

O PO CM O 

.-I iH 'T> CD 
vC th 



LIJ tx> Ol LIJ 
o rv tr% o 

CD »o o o 



ir\ 



CM 

.—I 

CO 



(Tk 



UJ 

CM 



r«v 

O 



Q 

UJ 



0 



Ui 



LTV Tk lf> 

CD CD O CD 
till 
L.U LU Ui UJ 
CM nO O 
o fNv t-i 

^ >0 CM >o 

nO tH tNO W 



O CM 



-x CD 



UJ 

If 

o 

o 

o- 



,H CM 
CD ^ 
•f I 

UJ UJ 

a ^ 

CD fv 
o CM 
,H CM 



PO iH O 
W tH .H CD 
1114 
UJ Ol UJ UJ 
M" O O CD 
CO .H .H O 
nD OO <X) CD t~ 

XIO.-I.HO CD 



oc 



cr 



o 



3 

u 

UJ 

X 

UJ 

M- *0 2- 

O 1—0 UJ 

I O I UJ 

UJ O UJ *13 

oo o o 

CM . 1 ^ in 

o O nO <t 

to 'O .H X 



UJ u>. in to 

(/) O CD CD CD 
CC I 1 • f 

UJ UJ U! UJ UJ 
Dv O nT -O CM 
Z fs. W IX CD 
— nO (M nO fx 
v4 PO tH 'O 
UJ - . «■ « 

X I I 



>— m .H CM fO 
t 4 O .H r 4 
— I I I 

Z UJ UJ UJ UJ 

UJ trv o .H 

O O CD CM CC 

— • o o ct> >o 

o w no m 

UJ • - - . 

X I 



•o 



t— 4 I I I 
Z UJ UJ UI UJ 

UI o On oo r-« 

O O r-l CM CM 

— . o eo rx .r 

O .H CM r4 

UJ • . - • 

X « 



in 

o 

o 



X 

3 



O 

o 



o 

i~ if\ m fx CO 

or o o o o 

<1111 
CL UJ UJ UJ UJ 
fX O CM < 
C*3 <■ fx o CM 

< .H 'O r- o 

X PO r4 nO CM 

1 I I 

UJ 
X 



o 

l~ -H CM CM 

X o r4 r4 t 4 

<^111 

X UJ UJ Ui UJ 

O u> IX vO 

_J o O PO fx 

< CD o w rx 

UI .H O r4 .H 

X . . . . 

UJ 



o 

I— W .H PO O 

X .H .H r4 O 

<1114 
X UJ UI UI UJ 

O ^ OD O 

CD tH .H <■ CD 

< CO CO O O 

X r 4 .H O O 

I I I 

UJ 

X 



UI 

o 



Z (Tk 



CM 

O 

I 

UJ 

o 

o 

O 

CM 



in 

o 

U- z 

o 

z 

z o 

X — 

3 I- 

3 CL 

o ro c 

C_) CD 

I z 

I- UI o 

Z fx — 

UJ <• 4- 



UJ 

o 



X 

OD 



o 

UJ 

o 

X 

o 



to 

a 

z 

o 

o 

UI 

to 

ITk 

ITk 

(3 

< 

to 



z 

X 

nO 




f 




X 

X 

\ 

CO 

h- 

7Z 

X 

6:^ 

<S? 

+ 

tu 

CP 

© 



LlJ 

X 

q: 

o 

CO 

t 

X 



Y-SCPLE - MBE-ai UMira'INCH. 

MIKLOSp TJ. JOB BOX M 

TRflhSIENT PLUS STEADY STATE 



4 

















CO 

ro 








LO 

CO 

CT 




1 

! 


Cv' 

1 ^1 

1 f! 


!- 

!cc> 

rO 




i 


i 


1 / 
'?o 

(O 

U<y 




* 




cs 

(S> 

G> ___ _ 




i 

j 

1 

! 


'i 

1 

1 


> 

CO 

CO 






c: 


CO 

\<s> 

■22 





6S0 



LU 



X o: 

Sh- 



CQ 



i: x: 
o o 

^ V 
CO CO 
I — I — 

zi3 =: 



CO 



CQ V 

°g 

LiJ 



CO 

CO 

Z] 



a 

! — i — 



^ LU 
CO ^ 

O CO 



a: 

q: 

m- 



<sl 

I 

U 



S03 



S00- 



ii03- 




X -SCALE = 1^0E+ea miTS/INCK 
Y-SCALE » 5.00E-01 miTS/IhICH. 

MIKLOSp TJ. JOB BOX M 

TRANSIENT PLUS STEADY STftTE 



I 

I 





1 


1 

CD 

CD 




i 




CO 

CD 

CD 






1 


i 




i 

1 


\ 1 


u 

\CD 






( 


1 




i 


\ 

( 


500 






\ 


CD 

'rD 


. 


i 

1 -- 




' 


r ' ' 

) 











S83 633^ zed- ?iS3- 



K-SCftLE - 1.83E+03 UNITS/INCH, 

Y-SCftLE «= 2.89E-01 UMITS/IMCH, 

MIKLOS, TJ„ JOB BOX M 

TRmSIEhT PLUS STEADY STATE 




I 




i 







thesM5816 

Analysis of a multi-degree of freedom vi 



3 2768 001 88333 3 

DUDLEY KNOX LIBRARY 










